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Abstract 

Based on magnetization, specific heat, magnetostriction, and neutron diffraction studies on 
single-crystal TbCo2B2C, it is found out that the paramagnetic properties, down to liquid nitrogen 
temperatures, are well described by a Curie- Weiss behavior of the Tb^"^ moments. Furthermore, 
below Tc= 6.3 K, the Tb-sublattice undergoes a ferromagnetic (FM) phase transition with the easy 
axis being along the (100) direction and, concomitantly, the unit cell undergoes a tetragonal-to- 
orthorhhombic distortion. For fields up to 90 kOe, no field-induced splitting of the Co 3d orbitals 
was observed; as such the internal field must be well below the critical value needed to polarize 
the Co 2>d subsystem. The manifestation of a FM state in TbCo2B2C is unique among all other 
isomorphous borocarbides, in particular TbNi2B2C (Tjv=15 K, incommensurate modulated mag- 
netic state) even though the Tb-ions in both isomorphs have almost the same crystalline electric 
field properties. The difference in the magnetic modes of these Tb-based isomorphs is attributed 
to a difference in their exchange couplings caused by a variation in their lattice parameters and in 
the position of their Fermi levels. 
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I. INTRODUCTION 



Whenever a family of compounds, containing 3d transition-metal (M) and rare-earth 
(R) atoms, manifests similar band structures, the trend in their magnetic properties can be 
rationalized in terms of the position of the Fermi level Ep within the density of states N{E) 
curve. The cubic, Laves-type RM2 family of compounds provides a best illustration.-'^'^ 
Here, it is commonly considered that the R-ion moments are localized while the M moments, 
if they exist, are itinerant. Assuming further that the mutual interactions among the mag- 
netic ions can be described within the molecular field theory and, in addition, using Landau 
phase-transition arguments together with Stoner and Wohlfarth-Rhodes criteria, then the 
magnetic properties of the whole RM2 family can be classified according to the first few 
energy derivatives of N{E) evaluated at Ep.—'^'^ This model is successful in explaining 
the evolution of the paramagnetic susceptibility, the type of the magnetic phase transition, 
the relative magnitude of the critical temperature, and, moreover, the magnetism of the 
3d subsystem; in particular it explains why the Ni sub lattice is magnetically silent in i?Ni2 
while the Co subsystem has an intermediate character: being strongly exchange-enhanced 
in RC02 {R = Y, Lu) while developing a delocalized moment in magnetic RCo2- For such 
RC02 series, it was found out that an adequate effective magnetic field at the Co site, H^£, 
is able, at and above a critical field H^^°, to induce a metamagnetic transition due to which 
the Co atom acquires an itinerant moment of ~l/iB- 

The above-mentioned interplay between the magnetism of M and R subsystems is ex- 
pected to be manifested also in the quaternary borocarbides RM2B2C family (see e.g. Ref. 
5 and references therein). Though the investigation of this interplay is of interest for the 
understanding of the magnetism (and superconductivity) of this family, there are only very 
few studies bearing on this subject and these are mostly limited to the i?Ni2B2C series. 
Nevertheless, the extensive magnetic characterization of this i?Ni2B2C series did clarify the 
influence of R (M ) on their overall magnetic (superconducting) properties;- in particular, it 
is estabhshed that while the M-doping hardly influences the magnetic properties (but leads 
invariably to a degradation of the superconductivity), the variation in R strongly modifies the 
magnetic structure: for R = Pr, Ho, Dy, the magnetic structures consist of a commensurate 
antiferromagnetic, AFM, state^S'^'^ that of i? = Tm is incommensurate modulated structure 
with 'g^~(0. 093, 0.093,0),-^'^ and those of i? = Er, Tb, Gd are incommensurate modulated 
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states with "g"^ (0.55,0,0).-'iiii^'i^'^ At Twf (which is <Tn), each of the latter states trans- 
forms into an equal-amphtude, squared-up state and in the particular cases of i? = Er and 
Tb, this transformation leads to a surge of a weak ferromagnetic component.— >^>^>^ It is 
an experimental fact that none of the Ni-based compounds manifests a ferromagnetic, FM, 
ground state though there is, at least for R = Pr, Ho, Dy, a strong FM intralayer coupling. 

It is recalled that the Ni subsystem in YNi2B2C is magnetically inactive^ while the Co 
subsystem in YC02B2C manifests an exchange enhanced paramagnetism:— based on the 
above model, these features indicate that Ep in the Ni-based compound is above a filled 3d 
band and, in addition, is not at a steep region of N{E) while for the Co-based compound the 
curvature of N{E) at Ep must be positive and nonnegligible. Indeed electronic structure 
calculationa^^i^'^i^ on i?Ni2B2C {R = Lu,Y) showed that Ep is situated at the top of a 
pronounced and narrow N{E) peak and this peak lies on the top edge of nearly filled Ni(3(i) 
bands: this, together with the intermediate electron-phonon coupling and the smaller Stoner 
factor, explains the surge of the superconductivity as well as the nonmagnetic character of 
the Ni-sublattice. Furthermore, in case of a FM order, such an N{E) peak at Ep would 
be exchange split due to the direct intra-atomic 4f-5d exchange interaction: such a split 
would be higher than the superconducting gap; it is an experimental fact that, due to their 
AFM structures, none of the magnetic i?Ni2B2C superconductors shows this splitting. 

Band structure calculations^^ on LUC02B2C showed that, for such isomorphous borocar- 
bides, the rigid band model yields a reasonable description of the band filling. Furthermore, 
Ep is situated at the decreasing but right-hand side of one of the peaks that receives a con- 
siderable contribution from the Co Srf-band, and that N{Ep) is of the same magnitude as 
that of LuNi2B2C: the latter finding is consistent with the observation that the Sommerfeld 
linear specific heat coefficients of YC02B2C and YNi2B2C are equal.— Furthermore, as that 
the first derivatives of N{E) are nonnegligible, then we expect the Co sublattice, similar to 
i?Co2, to develop an intermediate character (or even to be polarized) if H^^ > H^°: the 
intermediate character is indeed observed in YC02B2C;— the possibility of Co-subsystem 
polarization, on the other hand, would be addressed in this work. As far as the -R-sublattice 
magnetism is concerned, it is expected that the difference in the lattice parameters and in 
the electronic band structure of the Co- and Ni-based borocarbides would imply a modifi- 
cation in the character of the mediating RKKY-type interactions and consequently in the 
character of their magnetic ground states. 
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In this paper, we report on the extensive magnetic characterization of TbCo2B2C. The 
successful synthesis of a single-crystal sample of TbCo2B2C made it possible to identify 
unambiguously the paramagnetic as well as the ordered-state properties of the Tb sublattice: 
the former is dominated by a Curie- Weiss behavior while the latter is found to be a FM 
structure with the easy axis lying along the a direction. Such a FM state is in sharp 
contrast to the AFM-type mode of the isomorphous TbNi2B2C.— li^i^^iSiiSSiiSE This work also 
addressed the question of whether the surge of this FM state (with a strong Tb moment and 
consequently a strong H^^) is able to polarize the Co 3(i subsystem. Our results suggest 
that in spite of the intermediate character of the Co "id subsystem, the induced is not 
able to bring about an unambiguous spontaneous polarization. 

II. EXPERIMENT 

99.5% ^^B enriched polycrystals of TbCo2B2C were prepared by conventional arc-melt 
method. These polycrystals were used as feeding rods during a floating-zone synthesis, a 
process that we used for single crystals growth.— 

Characterization were carried out using magnetization [M(T, if), extraction method 
within the ranges 1.9 < T < 300 K and if < 90 kOe] and zero- field specific heat [C(T), 
semi-adiabatic method within the range 0.5 < T < 15 K, accuracy better than 4%]. A high- 
accuracy capacitance dilatometer— was used for measuring the thermal expansion or forced 
magnetostriction with a resolution better than 1 A. The relative change in length, measured 
along the cosine directions (/3i/52/33) when a field is applied along the cosine direction (xyz), 
is denoted as '^^^^^^A,,,(r) = [/(T, H) - kin, Ho))] /lo{To, Ho). 

Neutron-diffraction measurements were carried out at the Institut Laue-Langevin in 
Grenoble, France. Measurements were carried out on a powdered as-prepared arc-melt 
poly crystalline sample (since large quantity is desirable, see § III.D below) as well as on a 
single crystal sample. The powder diffraction patterns were collected within a temperature 
range 2 to 40 K using the DIB diffractometer with a selected incident wavelength of 2.42 
A; Rietveld refinements of both crystallographic and magnetic structures were carried out 
using the Fullprof package of Rodriguez-Carvajal (www.ill.fr/dif/Soft/fp). Single crystal 
diffraction studies were performed on the DIO four-circle diffractometer with A = 2.3606 A 
over a wide range of q space and within the temperature range 1.7 < T < 8K. 
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III. RESULTS 



A. Magnetization 




T(K) 

FIG. 1: (Color online) T-dependent Xdc (left-hand ordinate) and x'dc (right-hand ordinate) curves 
of TbCo2B2C measured along the a and c axes. The solid lines represent x = C/{T — 6); C 
and are the Curie- Weiss constants. The large symbols on the Xdc(^) curves represent the 
H/M('T, H ^ {)) as obtained from the Arrot plot. The upper-left inset shows, on expanded scales, 
the low-temperature XdciT, H/ /a) for H =1 kOe and 10 kOe. The lower-right inset shows the 
magnetization isotherms at T=1.9 K which demonstrates that for H >2 kOe, there is no hysteresis 
effect (see text). 

Figure [1] shows the T-dependent Xdc and Xdc curves that were measured at different fields 
and along the two principal tetragonal axes. The large symbols on the Xdc(^) curve represent 
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the H/M(T, H ^ 0) obtained from the Arrot plot: the excellent agreement emphasizes that, 
within this temperature range, the contribution of magnetic impurities is negligible. A Curie- 
Weiss fit down to liquid-nitrogen temperatures of XdciT, H\\c) gives yUcfr = 9.7(l)/iB and 9c= 
-29.4(1) K while that of Xdc{T,H\\a) gives /les = 9.6(l)/iB and ^^=14.4(1) K. Evidently, 
there are anisotropic forces but the effective moments are in excellent agreement with the 
value expected for a free Tb'^^ ion. Based on these anisotropic 6 values, the first Stevens 
coefficient in the crystal field description of a tetragonal symmetry is estimated to be i?2 = 
0.88(2) K; this compares well in sign and magnitude with that of TbNi2B2C (i?2 =1-2(1) 
K).— This similarity suggests that the crystalline electric field, CEF, at the Tb^^ site of 
both isomorphs are similar: indeed both sites have the same D^h symmetry and almost 
the same charge distribution. On lowering the temperatures toward liquid helium region, 
Xdc{T, H\\a) increases relatively fast and afterward tends toward saturation. Considering the 
characteristic magnetic features manifested in the magnetization, specific heat, and neutron 
diffraction (see below), this fast increase (which is followed by saturation) is caused by the 
process of approaching and the eventual onset of a FM order wherein the moments point 
along the a axis: Tc= 6.3(2) K is the point of maximum inclination. 

The upper-left inset of Fig. [T]shows, on an expanded scale, a magnetic hysteresis occurring 
at T < Tm =3.7(2) K and H=l kOe. The M{H\\a, 1.9 K) curve, shown in the lower-right 
inset of Fig. [T|, reveals that this hysteresis effect disappears for H > 2 kOe. Since this 
Tm-feature is sample-dependent (see below), it is attributed to a contaminating magnetic 
phase, the magnetization of which saturates completely to 0.05 yUe/formula unit for field 
higher than 2 kOe. Based on the weight ratios of the magnetic moments, the fraction of 
the Tb ions in this spurious phase relative to the major TbCo2B2C phase is estimated to be 
0.7%. 

Figure [2] confirms the above-mentioned magnetic anisotropy, due to which the a (c) is 
the easy (hard) axis. Within the studied ranges of H and T, the magnetization isotherms do 
not show any field-induced transition, rather, only a monotonic and steady increase (tending 
towards saturation) which is characteristic of a forced domain alignment: this supports 
the earlier inference of a FM order. Furthermore, the high-field magnetization increases 
as M{H\\a, 2K) = fi^p + XmH [/isp=7.6(3)/iB, Xhf=6.8(4) x 10~^ yUe/Oe corresponding to 
38(2) xlO~^ emu/mol] attaining /i(90 kOe, 2 K) = 8.2yUB- /^sp is only 3% lower than the 
reported moment of TbNi2B2C ( fiTh = 7.78 fiBr^) but 16% lower than the one expected 
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H (kOe) 

FIG. 2: (Color online) Magnetization isotherms of TbCo2B2C at different temperatures and for 
different field orientations. Two different set-ups were employed one with a field up to 90 kOe and 
another up to 50 kOe. The high-field magnetization evolves as M{H\\a, 2K) = -|- xmH and is 
represented by the solid line (see text). 



for a free Tb^"*" ion: this observed moment lowering confirms the above-mentioned influence 
of the CEF effects and that these effects are similar to the ones observed in the Ni-based 
isomorph. Alternatively, let us assume that the difference between /igp of TbCoaBaC and 
the reported moment of TbNi2B2C is due exclusively to the spontaneous polarization of the 
Co "id orbitals which are coupled ferrimagnetically to the FM Tb sublattice. Then, based on 
the relation /isp(TbCo2B2C) = /isp(TbNi2B2C) — 2./ico; the maximum possible Co moment 
would be /ico =0.25 /ie. Such a /ico value is surprisingly small, amounting to only 0.25 
hole/Co atom in the "id band; this value is almost one fourth of the Co moment encountered 
in the heavy members of the RC02 series,^ but it is similar to the value observed in the low 
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spin state of Co in, e.g., ErCo2.^^'^^ It is more likely that H^^ < H'^" (see § IV) since if there 
is any spontaneous Co polarization then the induced moment should be much higher than 
0.25 /iB- As the CW law describes well the paramagnetic susceptibility of TbCo2B2C (see 
above), then the exchange-enhancement factor for the Co-subsystem susceptibility must be 
extremely small. 

B. Specific lieat 

Figure [3] shows the zero-field magnetic specific heat and entropy of single-crystal 
TbCo2B2C obtained after subtracting the nuclear, electronic, and lattice contributions (the 
latter two were obtained from YCo2B2C).^^'^^ The nuclear contribution is of dominant im- 
portance only at very low temperatures and was evaluated from the diagonalization of the 
hyperfine Hamiltonian.— It is worth mentioning that at very low-temperature, the nuclear 
contribution is much stronger than the magnetic one: as such the propagation of errors 
due to successive subtraction of non-magnon contributions would eventually influence the 
absolute value of Cm(T); this may undermine the quality of the comparison between the 
theoretical and experimental magnon contributions (see Fig. E]). 

Both CmagiT) and SmagiT) curvcs do confirm the onset of the magnetic order at Tc. the 
former curve rises very sharply at while the latter manifests a pronounced change of slope. 
Considering the magnetic structure to be associated with the FM order of the Tb-sublattice 
(see below), we fit the experimental CmagiT) to the theoretical magnon expression of Eq. 
IA5I (see Appendix): as can be seen in Fig. [31 the excellent fit to Eq. I A5 1 gives the spin- wave 
stiffness coefficient i5=26.3(5) K and the gap parameter A=8.4(2) K. The high value of 
D is indicative of stronger effective exchange couplings. On the other hand, the value of 
A (which from Eq. IA3I is a measure of the anisotropic field) is consistent with the strong 
anisotropic features observed in the magnetization measurements. 

The inset of Fig. [3] compares the measured CmagiT) of a single-crystal sample with that 
of a polycrystalline one. Evidently the Tm-transition is sample- dependent: while CmagiT) of 
the polycrystalline sample manifests a pronounced event at Tni=3.6(2) K, that of the single- 
crystal hardly shows any anomaly. As mentioned above, the single-crystal sample contains 
only a 0.7% spurious phase and this concentration limit is lower than the resolution of the 
specific heat set-up. 
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FIG. 3: (Color online) Zero-field magnetic specific heat (symbol) and calculated magnetic entropy 
(dashed line) of single-crystal TbCo2B2C. The nuclear, electronic, and lattice contributions were 
already subtracted (the latter two were obtained from YC02B2C).— The solid line represents a fit 
to the expression of the magnon contribution (Eq. lASp assuming a FM order of the Tb sublattice. 
The inset (with the same scales as the main panel) compares the single-crystal magnetic specific 
heat (solid line) with that of a polycrystalline sample (symbol); the observed difference reflects a 
strong dependence on the sample form, contamination, and history. Nevertheless, their magnetic 
entropies approach each other for T > conflrming, as it should, the conservation of the total 
entropy (see text). 



C. magnetostriction 

Figure H] shows the forced magnetostriction isotherms measured along the a and h axis 
with H\\a. Once more (apart from the low-field, domain- wall sweeping-out region and the 
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saturated regime) there is no metamagnetic transition in these isotherms. The inset of Fig. H] 
indicates clearly that below T^, TbCo2B2C undergoes a spontaneous distortion which - based 
on Fig. m^a), see also § III.D - is attributed to an orthorhombic distortion of the tetragonal 
unit cell. A similar distortion was reported for TbNi2B2C.— '^'^'^'^ Then the behavior of 
the forced magnetostriction of FigJH can be understood as follows: at zero-field, there is an 
equal distribution of domains along each of the a and b axes; an applied field along, say, 
the a axis would involve a rearrangement of the orthorhombically-distorted domains and as 
such induces an increase in M{H\\a) [Fig. H] (b)], an increase in ^°°Aioo(-f^), and a decrease 
in °iOAioo(i^) [Fig. 11(a)]. 

D. Neutron Diffraction 

The thermal evolution of the powder diffractograms are shown Figs. [MSI For T > Tc, 
the patterns consist of the tetragonal crystal structure of TbCo2B2C and a small impurity 
phase. On the other hand, for Tm < T < Tc, it is evident that the magnetic reflections are 
piled up on the top of the nuclear Bragg peaks: a go=(000) mode. Considering that the 
paramagnetic state is dominated by the Tb moment, that the evolution of the isothermal 
magnetization and magnetostriction indicates no metamagnetic transition which can be 
related to the onset of Co moment, then this go mode must be related to the Tb sublattice. 
Furthermore, since the Tb ion occupies the special 2a site in the unit cell, then this mode 
must be FM: confirming the conclusions drawn from the magnetization, magnetostriction, 
and specific heat studies. Indeed, Fig. [U] (a) confirms that this pattern is a FM mode. 
Alternatively, if this mode is related to the Co sublattice then, due to the multiplicity of the 
4d site occupied by the Co atoms, the magnetic order should be either AFM (if only due 
to Co subsystem) or ferrimagnetic (if both subsystems are ordered): in the light of all the 
above-mentioned results, both possibilities must be ruled out (see also § IV). 

To extract more information, Rietveld analysis was carried out on the diffractograms 
measured within the range < T < Tc. Because of the structural distortion, we used the 
Immm space group together with the parameters given in Table [B Representative analyzed 
diffractograms are shown in Fig. [7] while the obtained cell parameters and the magnetic 
moment are given in Fig. [8l The analysis indicates that the lattice parameters [Fig. [8] 
(a-b)] undergoes a noticeable orthorhombic distortion below Tc which is consistent with 
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H (kOe) 

FIG. 4: Representative forced magnetostrictions (a) ^'''^Aioo(-f^, 3K) and ^^^Xoio{H, 3K) curves of 
TbCo2B2C are compared to (b) the isothermal magnetization at 3 K. The magnetostriction curves 
are given relative to their zero- field values. The inset shows the thermal evolution of the zero- field 
AL/L measured along the a axis. The arrow marks the Tc value which was determined from the 
specific heat measurement of Fig. 3. 

the magnetoelastic effects observed in Fig. HI Furthermore, the analysis revealed that the 
zero-field Tb magnetic moment is along the longest side of the base of the orthorhombic cell. 

Below Tin ~3.7 K, the thermal evolution of the diffractograms reveals two features [Figs. 
[6] (b and c)]: First, the intensity of the FM mode evolves smoothly and independently. 
Secondly, there is a surge of additional magnetic peaks [marked by the vertical arrows in 
the difference plots of Fig. [6] (b and c)]. Since Figs. [T]to|l]do not indicate any event that 
can be related to an order-to-order transition of Tb magnetic order, then these peaks can 
not be associated with the magnetic pattern of the Tb-sublattice. As these diffractograms 
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FIG. 5: (Color online) Representative neutron powder diffractograms of as-prepared polycrystalline 
TbCo2B2C. To aid in visualizing the magnetic modes, this plot is limited to temperatures below 
10 K and to scattering angle lower than 50°. The diffractogram at 6.4 K (denoted by thick solid 
line) can be taken as a demarcation between the paramagnetic and FM phases. 

were collected on as-prepared, arc-melt polycrystalline sample - which as mentioned above 
contains magnetic contamination - then these additional peaks are taken to be due to the 
same magnetic contamination, which is responsible for the hysteresis event in Fig. [T] and 
the weak specific heat anomaly in the inset of Fig. [S] (we argue, in § IV, that these events 
can not be due to an onset of a Co itinerant moment). 

Various g-scans within the range 1.7< T <8 K were performed on a single crystal of 
TbCo2B2C. A wide range of q space were scanned while maintaining the temperature con- 
stant at 1.7 K. Indeed, most of the nuclear (and ferromagnetic) peaks that satisfy the relation 
h + k + I = 2n were observed. In addition, we also looked, at 1.7 K, for any modulated 
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FIG. 6: (Color online) The magnetic diffractograms of as-prepared polycrystalline TbCo2B2C 
sample. These were obtained after subtracting either the pattern at 10 K (a and h) or the pattern 
at 4.2 K (c). (a) The FM mode; (6) the superposition of the FM and the magnetic contamination 
contribution at T=2.0 K (the strongest peaks of the later are denoted by vertical arrows); (c) the 
thermal evolution (T < Tm) of the magnetic patterns relative to the one at 4.2K. 

mode within a wide range covering —0.4 < h < 0.6, —0.4 < k < 0.4, and —0.3 < / < 0.3. 
Some weak reflections were observed. Fig. [9] shows the thermal evolution of one of these 
peaks and, for comparison, also that of the (0,0,4) peak. As evident, the intensity of the 
(0,0,4) peak, being due to magnetic and nuclear contributions, decrease smoothly and goes 
to the value of the nuclear intensity as the temperature reaches T^. In particular there is no 
visible variation in the (0,0,4) intensity when T is varied across Tm, indicating that the event 
at this temperature, does not belong to the main phase. On the other hand, the intensity 
of the other peak decays very fast as the temperature increases and is almost within the 
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TABLE I: Comparison of the space groups, atomic positions, and isotropic thermal parameters of 
TbCo2B2C and TbNi2B2C. The data on TbNi2B2C are taken from Lynn et al.^ The same thermal 
parameters reported for TbNi2B2C are also used for the analysis of TbCo2B2C. 

Tb Co B C 

TbNi2B2C PA/mmm position (000) {\^\) (0,0,0.357) {\\^) 
thermal factor 0.47 0.57 0.77 0.85 
TbCo2B2C Immm position (000) {\{)\) (0,0,0.354) (i^O) 

thermal factor 0.47 0.57 0.77 0.85 



experimental uncertainty when T > 2 K. Such a thermal evolution is similar to the features 
observed in the magnetization hysteresis (Fig. [T]), in the specific heat anomaly (inset of Fig. 
E]), and in the powder neutron diffractograms (Figs. [MZ])- This suggests that this, as well as 
the other weak peaks, are related to the same impurity phase as discussed above. Evidently, 
the ratio of the intensity of this contaminating peak to that of the (0,0,4) peak is extremely 
small in the single-crystal sample; in terms of the above mentioned impurity scenario, this 
means that the impurity concentration in the single- crystal sample is much smaller than 
the one in the polycrystalline case: this conclusion is supported by the observation that the 
intensity of the specific heat event at is hardly evident in the single-crystal sample (see 
the inset of Fig. [3]) . 

Based on the above conclusions, the powder diffractograms below Tm were analyzed, with 
the Rietveld method, as a superposition of three patterns: the nuclear, the FM mode, and 
a third unidentified contaminating magnetic phase. The Rietveld analysis of the nuclear 
and magnetic patterns of the main phase below is straightforward and gave a satisfacto- 
rily fit; the diffractograms are shown in Fig. [7] while the obtained structural and magnetic 
parameters are given in Fig. [HI the latter figure reveal that the lattice parameters evolve 
smoothly across Tm and, furthermore, the mismatch parameter {a — b)/a is the same as the 
one manifested in the magnetostriction experiment of Fig. H] (a). On the other hand. Fig. [8] 
(c) compares the reduced Tb magnetic moment, /i(T)//i(2 K), with the calculated Brillouin 
function, Bq{x). The total angular momentum quantum number is taken to be 6, repre- 
senting that of a free Tb^"*" ion. Since, as shown above, there are considerable CEF effects, 
then this Bq{x) curve should be taken as a lower bound; nonetheless, the overall thermal 
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FIG. 7: (Color online) Rietveld analysis on representative powder diffractograms. Within the range 
Tjn < T < Tc, the total contribution is composed of a nuclear and a FM component while for T < 
Tm it is a sum of three phases: a nuclear, a FM, and an impurity phase (see text). The thermal 
evolution of the lattice parameters and magnetic moment are given in Fig. [51 The Bragg ii-factor 
for the structural phase clusters around 2 while for the magnetic phase around 2.1. 



evolution of fi{T)/fi{2 K) follows reasonably well this 56(3;) curve; in particular, it reveals 
a smooth and monotonic evolution across the Tm region. Finally, the intensity contribution 
of the unidentified phase was calculated using the so-called profile matching^^ (or pattern 
decomposition) procedure. Since the crystal structure parameters of this unidentified phase 
are unknown, no significance should be attached to the fit of the impurity phase, only that 
all the additional weak peaks are associated with the impurity phase and that the presence 
of this phase would not modify the conclusions reached about the magnetic properties of 
the TbCo2B2C phase. 
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FIG. 8: (Color on line) The lattice parameters (a - b) and the reduced magnetic moment /i(T)/^(2 
K) versus the reduced temperature T/T^ (c). These parameters were obtained from the Rietveld 
analysis (see Fig. [7]). The solid line in the upper panel represent the Brillouin function for J=6. 

IV. DISCUSSION AND CONCLUSIONS 

The above results show that the paramagnetic susceptibilities are given by the CW be- 
havior of the Tb^"*" moments, that this contribution is strong enough to mask any magnetic 
contribution of Co-sublattice (if there is any), that /igfr and /isat are typical of Tb'^^ ion and 
are almost equal to the ones observed in TbNi2B2C (see Table [IT|) . and that the transition 
at Tc is identified as being due to the FM order of the Tb sublattice. In addition, the 
sample-dependent Tm-event is associated with magnetic contamination. Below we give a 
further argument in support of this latter identification. Obviously, if such a Tm-event is 
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FIG. 9: The upper panel shows, on a three dimensional plot, the thermal evolution of the (0,0,4) 
peak and an impurity peak. The latter together with its Gaussian fit is shown, on an expanded 
scale, in the inset. The lower panel shows the powder diffractogram at 2 K (taken from Fig. 7). 
The intensity of the (0,0,4) and impurity peaks (marked by a vertical arrow) are almost equal in 
the powder diffractrogram but are a factor of 20 different in the single-crystal diffractograms (no 
correction for multiplicity factors are considered). The wave length used for the measurement in 
the upper panel is 2.3606 A while for the lower panel is 2.42 A . 

intrinsic, then it must be either due to the Co- or Tb- sublattice. That the M{H,T < Tq) 
and ^^^Xioo{H,T < Tc) isotherms do not show any metamagnetic transition and that the 
same FM state of Tb sublattice is being maintained across Tm, then this Tm-event can not 
be related to a rearrangement (in direction or strength) of the Tb FM sublattice. 

Let us now discuss the claim that the Co orbitals which, being on the verge of ferro- 
magnetism, are spontaneously polarized. But the fulfilment of this possibility requires that 
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-f^eff — -^cr°5 ^ relation which can, apriori, be satisfied since H^g increases monotonically 
as T decreases [see Fig. [H](c)]. In this case, there should be two magnetic transitions: one 
related to the -R-subsystem and another to the Co-subsystem; just as in the case of, e.g., 
Ero.6Yo.4Co2.^'^' Then, as observed in other intermetallics,-'^'^'^ such induced Co-moments 
should couple ferrimagnetically with the Tb^"*" FM sublattice. As there are no spontaneous 
or field-induced metamagnetic transitions in the M{H,T < Tq) curves, then the possibility 
of a spontaneously polarized Co moment must be excluded; as such H^^ < H^j.'^{ see § III. A). 

TABLE II: A comparison of the lattice and magnetic parameters of the isomorphous TbCo2B2C 
and TbNi2B2C. The lattice parameters are reported for samples at LHe temperatures except the c 
parameter of TbNi2B2C; the later was estimated by normalizing its room temperature value using 
the values of HoNi2B2C..^ 
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Magnetic Mode 




TbCo2B2C 


3.535 


3.523 


10.560 9.7 6.3 


FM g=(000) 


7.6 
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3.554°^ 


3.534"^ 


10.44'' 9.8^= 15^ 


LSW, g=(0.45,0,0)^ 


7.78^*'^ 


Ref.-^, 






Ref..^, 


^ Ref.-^ 





It is significant that the FM structure of TbCo2B2C is drastically different from any 
of the reported magnetic structures of borocarbides,- in particular TbNi2B2C even though 
these Tb-based isomorphs are similar in most (if not all) of the single-ion CEF-influenced 
properties such as the anisotropy, the strength, and the orientation of Tb moments (see § 
III. A and D). However, these isomorphs are distinctly different (see Table HTj) in the value 
of their transition points, in their magnetic structures, and in the overall features of their 
H — T phase diagrams (a simple one-boundary FM phase versus a cascade of field-induced 
phase transitions). These differences suggest that the effective Tb-Tb magnetic couplings 
must be different and, furthermore, those nesting features^^ that are responsible for the 
modulated mode in TbNi2B2C must be absent in TbCo2B2C. These differences, prompted 
by the introduction of the Co atoms, suggest that the configuration of their electronic 
structures [in particular the position of Ep within the N{E) curve and the generalized 
susceptibilities] must be different. As mentioned in § I, these arguments are consistent with 
the findings of the electronic structure calculation on LUC02B2C:— these calculations provide 
an explanation for the surge of the enhanced paramagnetic character of YC02B2C and the 
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absence of superconductivity in any of the i?Co2B2C compounds (even though YC02B2C 
has the same Debye and Sommerfeld coefficients as YNi2B2C).— 

Finally, it is recalled that the indirect exchange coupling in metallic magnets are usually 
written as:— 



F{x) 



X cos(x) — sin(x) 



{xr 

where F is the s — f exchange coupling, n is the carrier concentration [governed by by 
N{Ef)], Rni is the distance separating the moments, Ep and kp are the Fermi energy and 
wave vector, and A {> R) is the mean free path. Considering that these couplings manifest 
a quadratic dependence on N{Ef) and a sensitive sinusoidal dependence on the moments 
separating distances, then it is no surprise that the combination of difference in the electronic 
structure^^ and in the lattice parameters (see Table [T]) would lead to strong variation in the 
magnitude and sign of the coupling constants and as such to drastic difference in the magnetic 
structures of these Tb-based isomorphs. In fact this difference is not limited to these Tb- 
based isomorphs, our preliminary studies on the magnetic structures of the -RC02B2C series 
showed that this is valid for the whole -RC02B2C magnets:— as an example, the FM mode 
is observed in TmCo2B2C and H0C2B2C. 
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APPENDIX A: MAGNON CONTRIBUTION FROM THE FERROMAGNETIC 
TB-SUBLATTICE OF TBCO2B2C 

The above-mentioned ferromagnetic order of the Tb-sublattice suggests that the dominant 
exchange interactions within the same layer (approximated by a positive Ji) as well as those 
among different layers (approximated by a positive J2) are ferromagnetic. Let us assume 
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that the main contributions to the magnetic energy is due to the above-mentioned exchange 
couphngs and anisotropic crystalhne electric field interactions. Within the low temperature 
regime of interest, the anisotropic interactions can be approximated by an effective field Ha 
(in energy units) which, for this particular case, forces the moments to points along the a 
axis. The Hamiltonian (under zero external field) can be written as: 

n = - Yl -J^^^- - E -J^^^ ■ -Ha.J^s", (Al) 

i,jeA,B <ij>,ieA,jeB ieA,B 

All symbols have their usual meanings. The first term sums the bilinear products of two 
neighboring spins of the the same layer A and afterwards the contribution of all layers 
are added together. The second term sums all the bilinear product of two neighboring 
spins (each belong to a different but an adjacent layer). The third term is a sum over all 
single-ion anisotropic energies. Using standard linear spin-wave approximation (considering 
a non- interacting magnon gas), Eq. lAll can be diagonalize to give the following dispersion 
relation: 

= SJi{A — 2 cos bk^ — 2 cos akz) + 



cos{r,kz ~^ 9^1/) cos(r)kz ~\~ nkx r,k^ 



ai 



,2"'z 2 X ' 2 y> ^^^\2"'Z I 2 X 2 V> 

where z (the quantization axis), x, y axes are, respectively, along the a (easy direction), 6, c 
directions of the crystallographic unit cell. In the long wave limit, up to second order, this 
simplifies to: 

fitufc = A + Cxkl + Cykl + Czkl. (A2) 
Assuming a weaker orthorhombic distortion (a~b and Cx ~ c^), one gets: 
Cx = Sa^^Ji + J2) = Cz 

Cy = SC^J2- 

The energy gap {k = 0) is: 

A = Ha (A3) 

The expression for the density of states (obtained from integrating over the constant 
energy surface, e = u) is: 



~ (27r)3 J I Vhujk I " 27r2 Cx^ya^c ^ ^ 
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The magnon contribution to the total energy is: 



where Eq is a constant independent of temperature. Using Eqs. IA2[ IA31 and IA41 the molar 
specific heat is {u is the number of moles): 



IdE 



R 



V 



where 



The final result can be rearranged to give: 



Cmag{T) — 3 ^„ 

8.7r2ij2 



15.i?Ai ^ , nA, 



n=l 



T 



15V^j ^V^j ^U^J 



(A5) 



where 



D = 23 5(Ji + J2)«^2' 



(A6) 



is the spin-wave stiffness coefficient^ which is a measure of the effective coupling strength. 
This expression reproduces the high temperature limit (T >> A), 

CM(T)oce~^/^T3/2, 

which for A = gives the well known |-Bloch expression.— At lower temperatures (T << 
A), the strong dependence on the gap is emphasized by the expression: 



Cm{T) oc 



23 



